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Abstract 



We present an update of the One-Meson-Exchange (OME) results for 
the weak decay of s- and p-shell hypernuclei Q, paying special attention 
to the role played by final state interactions between the emitted nucleons. 
The present study also corrects for a mistake in the inclusion of the K and 
K* exchange mechanisms, which substantially increases the ratio of neutron- 
induced to proton-induced transitions, Tn/Tp. With the most up-to-date 
model ingredients, we find that the OME approach is able to describe very 
satisfactorily most of the measured observables, including the ratio Tn/Tp. 

PACS number(s): 21.80.+a, 24.80. +y, 13.30.Eg, 13.75.Ev, 13.75.Cs 

I. INTRODUCTION 

Hypernuclei are bound systems of non-strange and strange baryons. In present facilities, 
hypernuclei are created with hadronic reactions -as (vr,K) at Brookhaven and KEK- or 
electroproduction ones -as (e,e'K) at TJNAF-. The decay of those objects proceeds via the 
weak interaction which is 10 orders of magnitude slower than the strong one and violates 
parity, isospin and strangeness. For the very light hypernuclei the mesonic decay mode 
(A vrA^) is dominant, but as A increases so does the Pauli blocking acting on the outgoing 
nucleon, and hypernuclei mainly decay via the one-nucleon induced non-mesonic mode, 
AA^ — i> A^A^. Since the pioneering phenomenological model of Block and Dalitz 0, many 
approaches have been developed to understand the dynamics of the decay, and the results 
have been collected in extensive review articles &-§[. Many works are based on a meson 
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exchange picture, either using a simple one-pion-exchange mechanism |^J^, or including also 
heavier meson exchanges, such as the p or the complete pseudoscalar and vector meson 
octets The effect of correlated-two pion exchange mechanism [p!0|-[T^ as well as the 

role of A/ = 3/2 transitions, implemented in a meson-exchange picture |T^, have also been 
investigated. A four quark weak transition effective Hamiltonian, corrected by QCD, which 
contains both A/ = 1/2 and 3/2 transitions, has also been applied in the study of the weak 
decay of hypernuclei [^^-|T^. In general, the weak decay rates are reasonably reproduced 
by the models but the ratio of neutron induced {An — > nn) to proton induced {Ap — > np) 
decays, r^/Fp, turns out to be smaller than the experimental value, which is of the order 
of 1 or larger |T^-^, although new recent theoretical progress has been achieved into the 
solution of this puzzle |ll2|,[T6|. In addition to AA^ NN, the decay can also proceed via the 



two-nucleon induced process, ANN NNN, originally studied in Ref. [^, which amounts 
to about 15% of the total width p2|-|2^ and cannot be neglected in the experimental analysis 
trying to extract the ratio Tn/Tp. 

The high momentum of about 400 MeV/c transferred in the AA^ — > A^A^ reaction makes 
this process quite sensitive to short range physics. As a consequence, the strong baryon- 
baryon interaction both in the initial and final states plays a quite important role. The 
purpose of the present work is to revisit the one-meson-exchange model of Ref. , with an 
especial interest in quantifying the effect of final state interactions, as well as the uncertain- 
ties of different model ingredients. In doing so, we will also point out a sign error encountered 
in certain transitions mediated by the exchange of strange mesons which, when corrected, 
gives rise to a considerably increase in the Tn/Tp ratio. We will analyze the sensitivity of our 
results to the way final state interactions of the emitted nucleons are implemented, using 
different prescriptions to obtain the A^A^ scattered wave function. Our study concludes that, 
with the appropriate treatment of final state interactions and with the correct sign for the 
contribution of the strange mesons, the one-meson-exchange model is able to describe very 
satisfactorily most of the measured observables, including the elusive ratio Tn/Tp. 



II. WEAK TRANSITION POTENTIAL 

The weak transition potential is obtained by following the model of Ref. [|l| . In analogy to 
One-Boson-Exchange (OBE) based models of the strong interaction, the present formalism 
includes not only the exchange of the long-ranged pion, but also more massive mesons 
which account for shorter distances. This potential has been presented in previous papers 
and, therefore, it is not going to be discussed here in a great detail. However, we would like 
to show its expression in coordinate space, which represents a compact way of including all 
the mesons and transition channels in the mechanism. The r'-space potential then reads: 

= EEKl^nr) = EEKf^nr)Oa(r)/« , (1) 

i a i a 

where the index i runs over the different mesons exchanged {i = 1,...,6 represents 
7i,r],K,p,uj,K*) and a over the different spin operators: central spin independent (C), cen- 
tral spin dependent (SS), tensor (T) and parity violating (PV). The angular dependence is 
represented by the Oa{f) operator, explicitly given by: 
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OJr) 



i C (only for vector mesons) , 

c?icr2 SS , 

Sl2{r) = 3(Tif(T2f - (71(72 T , 

i (72r PV (for pseudoscalar mesons) 

[(Ji X (T2] r PV (for vector mesons) , 



while the isospin operator, J*^*\ takes the form fi T2 for isovector mesons (iTjp), 1 for isoscalar 
mesons {r],u!) and a combination of both operators for the isodoublet {K,K*). 

The different pieces of the potential, V^^\r), are found by Fourier transforming the 
potential in momentum space. In particular, for pseudoscalar mesons the g-space potential 
reads: 
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where Gprn^ = 2.21 x 10 is the Fermi weak constant times the pion mass squared, 5'bbm(») 
the coupling at the strong baryon-baryon-meson (BBM^*)) vertex, A^*) and 5^^) the parity 
violating (PV) and parity conserving (PC) weak couplings respectively, yU*-*^ the meson mass, 
and M2 (Ml) the average of the baryon masses at the strong (weak) vertex. 

We want to note here that the convention has always been to direct the momentum 
towards the strong vertex and, in connection to this, care must be taken when this expression 
is applied to K {K*) exchange, since in this case vertex 1 refers to NNK{K*) and vertex 
2 to ANK{K*). Therefore, the combination of the strange meson transition amplitudes 
with the other non-strange meson contributions requires to exchange, in the former ones, 
the space, spin and isospin quantum numbers between the A and N in the initial state 
and between the two emitted nucleons in the final state. For a particular transition from an 
initial AN{Lq, Sq, Tq) channel to a final NN{L, S, T) one, this exchange operation introduces 
the additional factor: 
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where ti = t2 = = ^ is the nucleon isospin and the isospin assigned to the A isospurion, 
which can take the values | or | depending on whether one is considering a change in isospin 
of 1/2 or 3/2. The way to implement A J = 3/2 transitions in the present mechanism can 
be found in Ref. |]13[. The factor in Eq. (|^) is nothing but a sign in front of the PC and/or 



PV amplitudes, which is: 
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It is common to write the total decay as the sum of definite spin-isospin amplitudes. For 
example, for s-shell hypernuclei one has B: 
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If we adopt the A/ = ^ rule - which we always do for pseudoscalar meson exchange 
|T3| , the kaon and pion potentials differ in a minus sign only in the / amplitude in the 



only PV amplitude with S = So{= 1). 

In order to account for finite size effects we include a monopole form factor at each 
vertex, = (A^*) — /x*^*) + g^), where the value of the cut-off, A^, depends 

on the meson (/xW). The updated expression for the regularized potential was given in Ref. 



25 1 . In previous calculations we used the cut-offs given by the Jiilich B interaction [26 



The reason was that unlike the early YN Nijmegen model [|2^, which used different cut 



offs depending on the irreducible representation of the baryon-baryon channel, the Jiilich B 
was the only OBE model which used different cut-offs depending on the meson. However, 
the Nijmegen group has recently made available new baryon-baryon interactions in the 
strangeness S = 0, — 1, — 2, — 3 and —4 sectors pH|. These potentials are based on SU(3) 



extensions of the models in the S = and —1 sectors, which are fitted to the experimental 
data. The authors of Ref. give six different models, which fit the available A^A^ and YN 



scattering data equally well but are characterized by different values of the magnetic vector 
F/{F + D) ratio, ranging from 0.4447 (model NSC97a) to 0.3647 (model NSC97f). The 
advantage of these new models is that the form factors depend on the SU(3) type of meson. 
The momentum space potential for each meson is multiplied by the regularizing factor of 
Gaussian type, exp(— g^/A^), with a cut-off Ai for the singlet meson, Ag for the non-strange 
members of the meson octet and A;^ for the strange meson. In order to accommodate to 
our own formalism, which uses a monopole form factor at each vertex, we will match the 
Gaussian to a function of the type [A^/(A^ + g^)]^ at |g| ~ 400 MeV/c, the most relevant 
momentum transfer in the AA^ — > A^A^ process. Since, by definition, both functional forms 
give 1 at g* = 0, our alternative expression with the modified cut-offs A, listed in Table |, 
gives an excellent reproduction of the Gaussian NSC97 form factors up to a momentum 
transfer of about 600 MeV/c. 

As it is well known, one of the sources of uncertainty in OBE models comes from the 
coupling constants between baryons and mesons. In the strong sector the different inter- 
action models use SU(3) in order to obtain the BBM^*-* couplings that are not constrained 
experimentally. In the weak sector, only the decay of the A and S hyperons into nucle- 
ons and pions can be experimentally observed. For the other mesons, SUyj{Q) represents a 
convenient tool to obtain the PV amplitudes, while for the PC ones, we use a pole model 
with only baryon pole resonances. Details of how these coupling constants are derived 
can be found in Ref. H and, in the recent work on the weak decay of double-A hypernuclei 
29| , one can find an updated list of the S-wave (PV) and P-wave (PC) coupling constants. 



the latter obtained from the pole model using different parametrizations of the strong BB 
interaction. 
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A measure of the amount of parity- violation in the weak decay, is given by the asymmetry 
in the angular distribution of protons coming from the decay of polarized hypernuclei. This 
asymmetry is given by: 

^ 'J+1 Tr{MM^) 

where Py is the polarization of the hypernucleus, characteristic of the production reaction, 
and Ap the hypernuclear asymmetry parameter, characteristic of the weak decay. In the ex- 
pression above, J is the spin of the hypernucleus, Sy is the J-spin operator along the direction 
perpendicular to the reaction plane, and Ad the hypernuclear transition amplitude. In order 
to compare with experiments, one has to multiply Ap by the model dependent quantity Py, 
which has to be determined theoretically for each hypernucleus By using a shell-model 
for the initial hypernucleus and assuming spherical configuration with no mixing, one can 
express this amplitude in terms of two-body transitions, AA^ — >■ A^A^. The dependence of 
weak decay observables on the deformation of the initial (p-shell) hypernucleus was investi- 
gated in Ref. 0], by means of the Nilsson model with angular momentum projection. It was 



found that deformation effects change these observables by about 10% from the spherical 
limit, a deviation that although non-negligible it is smaller than the present experimental 
uncertainties. 



III. EFFECTS OF THE STRONG INTERACTION 

Because of the lack of stable hyperon beams, access to the AS = — 1 AA^ interaction 
is limited right now to the decay of hypernuclei. Hence, extracting information of the 
elementary weak two-body interaction requires a careful investigation of the many-body 
nuclear effects present in the hypernucleus. 

On the one hand, one must consider that the interacting nucleon and A hyperon are 
bound in the nucleus and they should be described by bound-state single-particle wave 
functions, obtained from appropriate mean-field or Hartree-Fock potentials. Note, however, 
that since the mass excess of 176 MeV in the initial state is converted into kinetic energy of 
the final particles, the nucleons emerge with a large momentum of about 400 MeV/c and the 
decay is not very sensitive to the details of the single-particle wave function. In Ref. it is 
shown that the decay rates obtained from various realistic A wave functions differ by at most 
15%. On the other hand, the large momentum transfer implies that the AA^ NN decay 
process is very sensitive to the short-range correlations induced by the strong interaction. In 
the initial system, one must then replace the mean-field two-particle AA^ wave function by a 
correlated one that accounts for the effects of the strong YN interaction at short distances, 
which are not considered in mean-field models. Correlated wave functions, obtained from 
the soft and hard core Nijmegen YN interactions |]27| , |32[| by solving the corresponding finite 
nucleus scattering amplitude (G-matrix) were compared in Ref. The differences 
between the wave functions obtained with the two models are already significant below 
0.75 fm and give rise to decay rates that differ by slightly more than 15%. However, since 
the present YN models are not constrained enough by the scattering data to resolve this 
discrepancy, we will admit this uncertainty in our initial AA^ wave function. This uncertainty 
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also justifies the use of a spin-independent parametrization for the AA^ correlation function 
which, when multiplied with the uncorrelated one, was shown to give a decay rate in between 
those using the soft and hard core correlated wave functions . 

Any realistic calculation must also take into account the fact that the two nucleons 
emerging from the decay feel their mutual influence, as well as that from the residual (A- 
2)-particle system. However, as mentioned before, the most important contribution to the 
decay comes from the so-called bacJc-to-bacic kinematical situation, in which the two nucleons 
emerge with the largest possible momentum of around 400 MeV/c. For these fast moving 
nucleons the distortions with the residual nucleus should be small, and the importance of 
such effects are further diminished for inclusive observables such as the decay rates. In other 
words, one could take into account the interaction of the emitted nucleons with the residual 
nucleus through an optical potential. However, the real part should play a minor role at 
those high energies, and the imaginary part will remove flux from the A^A^ channel which 
will reappear in other multi-nucleon channels, such that the total strength or decay rate, 
which is the quantity we are interested in, will not be modified. We would face a completely 
different situation if we were interested in calculating the energy distribution of the nucleons, 
in which case one should consider these final state interaction effects, which can change the 
energy, direction and charge of the primary nucleons emitted in the weak AA^ — > A^A^ 
transition, producing as well low energy secondary nucleons. The Monte Carlo simulation 
of Ref. p3| finds indeed that the final state interactions affect the nucleon distributions 



mostly at energies below 50 MeV. In the particular case of \^C, it is found that final state 
interactions produce roughly 1/3 more protons at those low energies, half of which are due 
to charge-exchange reactions, but barely affect the distribution at higher energies and leave 
the decay rate intact. 

Since in the present work we will only be interested in the decay rates, we will limit the 
treatment of final state interactions to those related to the mutual influence between the two 
emitted nucleons. In contrast to the strong AA^ interaction, the nucleon-nucleon one is much 
better constrained by the huge amount of scattering data and, although differences in the 
wave function can also be observed with various potential models, they are only significant 
below 0.5 fm, having a moderate influence in the decay rates, as we will see. 

In the literature, one finds a variety of ways of dealing with the FSI between the emitted 
nucleons. Some works do not include FSI, others use a phenomenological correlation function 
that simply multiplies the uncorrelated one, and others use various approximations to the 
A^A^ scattering equation. In the next section we will analyze the sensitivity of the non- 
mesonic decay observables to these various choices of FSI. The details on how the scattering 
A^A^ wave function is obtained from the Lippmann-Schwinger (T-matrix) equation, as well 
as from other simplified approximations, are given in the appendix. 



IV. RESULTS AND DISCUSSION 

We start this section by showing the weak decay observables for ^He and obtained 
with the model parameters used in Ref. [Q, but correcting the sign in the appropriate K 
and K* amplitudes, as pointed out in Sect. II. Hence, we use the coupling constants of 



the NSC89 strong YN interaction together with a monopole form factor of the type 
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F{q'^) = (A^ — yU^) / (A^ + g*^) at each vertex. Each meson has a different cut-off value, A, and, 
since this particular Nijmegen model distinguishes cut-offs only in terms of the transition 
channel, we choose, as in Ref. [|1|, the values of the Jiilich YN interaction |]26|. The results 
are collected in Table |I| and we remind here the main results given in Ref. [|l| for ]^C, namely 
a non-mesonic rate of Tnm = 0.75 Fa, where Fa = 3.8 x 10^ s~^ is the decay rate of a free A, 
and a neutron-to-proton ratio F^/Fp = 0.068. 

As we see, the corrected model gives a slightly larger decay rate, together with a remark- 
able increase in the neutron-to-proton ratio. This is due to the fact that the corrected sign 
gives rise to a constructive interference between the vr and K contributions to the 1=1 
^Si Pi PV amplitude, instead of the destructive one found in Ref. |Q and displayed in 
Fig. 8 of that work. On the other hand, the incorporation of K exchange to the vr mechanism 
lowers the PC amplitudes. This was pointed out in Fig. 7 of Ref. and Figs. 4 and 5 of 
Ref. ||3^, where the contributions of all different mesons are displayed, the vr and K ones 
showing a destructive interference in the PC amplitudes. The present constructive interfer- 
ence in the PV channels is especially relevant in the neutron-induced rate, F„, which lacks 
from the PC tensor transition for L = initial states. Conversely, the proton-induced rate, 
Fp, is dominated by the PC tensor amplitude, hence its value is lowered when K exchange 
is added to the vr mechanism, although the reduction will now be more moderate due to 
the constructive interference in the less dominant PV transition. Altogether, the correct 
incorporation of K exchange produces, in comparison to the results given in Ref. a more 
moderate decrease in the total non-mesonic decay rate, hence a slightly larger rate is now 
obtained, and a Fn/Fp ratio which is almost a factor of 3 larger. We note that these effects 
of K exchange have also been pointed out in Ref. |]16[, where the vr and K mechanisms are 
considered together with a description of the decay in terms of quark degrees of freedom, 
and in the recent work of Ref. |T2[, where, in addition of vr and K exchange, the role of 
correlated two-pion exchange was also studied. 

Results for the asymmetry of the emitted protons from the weak decay of polarized 
hypernuclei are also listed in Table Since this is an observable tied to the interference 
between PV and PC amplitudes, it will also be influenced by the change of sign in the 
abovementioned PV K and K* transitions. Indeed, the magnitude of the asymmetry Ap 

if 



has increased almost a factor of 2 from the values given in Ref. [H which are —0.27 for \B.e 



and 0.16 for fC. 

The fact that the new Nijmegen potentials include different form-factors depending 
on the meson, makes it possible to treat all the strong interaction ingredients involved in 
the weak decay transition AA^ NN consistently within the same model. We remind that 
the strong interaction plays a role not only in the strong vertices and form factors, but also 
in the PC piece of the weak vertex through the pole model and in the determination of the 
distorted final state NN wave function. In Table |T| we show the decay observables for ^He 
and obtained with the models NSC97a and NSC97f. 

We observe substantial differences between both models. These might come from the 
different coupling constants used in the strong and weak vertices of the transition AA^ — > A^A^ 
amplitude, or from the different distorted wave functions. To disentangle both effects, we 
show in Table |V| the results for \B.e, obtained using the same strong interaction model 
for the vertices, which we choose to be the NSC97f one, and different distorted NN wave 
functions obtained with four NN interaction models, namely the NSC97a and f EH], the 
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NSC93 and the Bonn B potential The resulting NN wave functions in various 

relevant channels are shown in Fig. |l| for a relative momentum = 386 MeV/c. Some 
minor differences appear only in the short distances, which should affect the decay rates 
moderately. By inspecting the results of Table one can indeed see that the non-mesonic 
decay rates in ^He vary from 0.3 to 0.4 Fa, the neutron-to-proton ratio F„/Fp from 0.45 
to 0.5, and the asymmetry from —0.7 to —0.6. In conclusion, it is mainly the differences 
associated to the strong couplings and their influence on the weak ones the origin for the 
differences observed in Table |T|. Going back to this table, we see that there are already some 
differences between the results obtained with models NSC97a and NSC97f when the K meson 
contribution is added to the vr meson one, but this cannot explain the differences in the final 
results when all mesons are included. Compare, for instance, the nonmesonic decay rate in 
^He of 0.43 obtained with the NSC97a model versus the value 0.32 obtained with the NSC97f 
one. In fact, by analyzing the weak couplings listed in Ref. [^, we can see some appreciable 
differences for the NNK coupling constants, which are larger in the NSC97f model, thus 
enhancing the destructive interference of the K meson and producing a reduced rate in this 
case. We also observe that the weak vector PC NNuj coupling constant obtained from the 
NSC97a model is one order of magnitude larger than that obtained from the NSC97f one. 
This is a consequence of a subtle cancellation in the pole model expression, which involves 
the difference between the strong NNuj and AAo; coupling constants. Although both models 
use the same strong NNuj couplings, the differences in the AAu ones induced by the freedom 
allowed in the magnetic F/ {F + D) ratio may give rise, as is the case here, to substantial 
differences in the values of the PC weak coupling constants. The larger value found in the 
case of the NSC97a model enhances the constructive interference of the uj meson contribution 
and explains the final difference between the decay rate obtained with the two Nijmegen 
models. This discussion makes also clear that one must admit this amount of uncertainty on 
meson-exchange models describing the nonmesonic weak decay of hypernuclei. We should 
therefore conclude that the nonmesonic decay rate for ranges between 0.55 and 0.73 Fa, 
the neutron to proton ratio F„/Fp between 0.29 and 0.34, and the hypernuclear asymmetry, 
more stable, is 0.36. In AHe, the rate obtained lies in the range 0.32 — 0.43 Fa, F„/Fp in the 
range 0.34 — 0.46 and the hypernuclear asymmetry parameter is —0.68. We observe that our 
results for the nonmesonic decay rates compare satisfactorily with the experimental data, 
especially after considering that the size of the nonmesonic two-nucleon induced mechanism, 
not included here, is around 0.3 Fa for medium to heavy hypernuclei ||2^j2^ . The ratio F„/Fp 
is considerably larger than that found in our previous works and it now lies practically 
within the lower side of the experimental errors. We note that a larger ratio could have been 
obtained, as we show below, if we had used a phenomenological treatment of FSI effects or 
we had ignored them althogether. However, we want to stress here that the results shown 
in Table JTl| have been obtained consistently within the same strong interaction model and 
that involves also the dynamical generation of the corresponding A^A^ wave function through 
a Lippmann-Schwinger equation. 

The asymmetry parameter Ap has increased by more than a factor of two with respect to 
our previous works. As it has already been mentioned in connection to the results of Table ||, 
this is a direct consequence of the change of sign in the K and K* exchange PV amplitudes 
for transitions that do not change the spin S, such as / in Eq. (H). Indeed, this amplitude 
interferes with c and d and the magnitude of the latter tensor amplitude is of especial 
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relevance in meson-exchange models, as the one employed here. The increased asymmetry 
values compare less favourably with the experimental data. While the result for is 
still within errors, the large negative value obtained here for ^He is very far away from the 
positive small result obtained in Ref. We note that, using a weak coupling scheme, one 
can relate the hypernuclear asymmetry parameter Ap to that associated to the elementary 
reaction AA^ — > A^A^, a\ (see Eq. (9) in Ref. [Q]). In the case of ^He, oa = Ap and for 
the relation is oa = —'2Ap. Therefore the parameter qa of the elementary reaction that can 
be extracted from our results in Table turns out to be very similar in both hypernuclei. 
However, the present experimental results seem to contradict this finding because a small 



positive value is obtained for ^He p6[, while a large negative one was found in [38 



In Ref. 1 36] arguments are given to try to understand these differences on the basis of the 
fact that ^He is a s-shell nucleus and a^C a p-shell one, the latter one allowing transitions 
from initial AA^ relative P-wave sates. However, it was already noticed that most of the 
decay for p-shell nuclei comes from relative S'-wave initial states and we have checked 
here that removing the initial P-wave amplitudes, the decay rate of using the NSC97f 
model gets slightly reduced to 0.48 and the asymmetry parameter Ap only changes to 0.34. 
Hence, meson-exchange models do not explain the present experimental differences between 
the elementary asymmetry for ^He and ]^C, and this seems to be a new puzzle in the study 
of the weak decay of hypernuclei. This poses a challenge to the theoretical models and calls 
for further experimental efforts. Besides valuable new information from hypernuclear weak 
decay experiments, it is also possible to gain direct information from the inverse reaction 
pn pA, already studied theoretically [ p5| , ^^ -^ and presently under preparation at RCNP 
(Osaka). 

One of the purposes of this paper was to study the influence of FSI in the AA^ NN 
transition. The results of Table |I^ have already shown that, using NN wave functions 
obtained from various realistic NN potentials, one obtains similar values for the weak decay 
observables. However, one finds in the literature several possibilitites for the treatment of 
FSI. They are either ignored, or included through a phenomenological function, such as 
/phen = 1 — joiner) with = 3.93 fm~^, multiplying the free NN plane wave, or through 
a distorted wave function built up from the i^'-matrix rather than from the T-matrix. In 
Fig. we show the resulting NN wave functions in various channels for a relative momentum 
Pr = 386 MeV/c, obtained from the various approaches to FSI. The solid lines represent 
the real and imaginary parts of the wave function obtained from the proper solution of the 
Lippmann-Schwinger equation (T-matrix). The imaginary part is the minoritary component 
in all panels. The dotted line represents the result obtained when only the real part of the T- 
matrix is retained in the construction of the distorted wave function from Eq. ( |A.11| ), which 



is further approximated by its principal part value. This was the prescription followed in 
Ref. [0. The dashed line corresponds to a A^-matrix approximation, which would be the 
appropriate solution for standing waves, i.e. non-propagating solutions, as is the case in the 
nuclear medium, where Pauli blocking prevents two interacting nucleons of the Fermi sea 
from exciting intermediate NN states that are already occupied by other nucleons. Finally, 
the dot-dashed line displays the product of the phenomenological correlation function with 
the corresponding non-interacting Bessel function. We observe important differences beween 
the various approaches to the treatment of FSI, especially in the case of a phenomenological 
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correlation function. Obviously, this will have an influence on the decay observables. As an 
example, we show in Table ^ the results for ^He, obtained using the NSC97f model and 
taking all these different treatments for the distorted NN wave functions. The rates can 
differ substantially, by more than a factor of two in some cases. The neutron-to-proton ratio 
can also be flctitiously enhanced by using the phenomenological approach to FSI effects, or 
no FSI at all, while the asymmetry remains more stable. 



V. CONCLUSIONS 

We have revisited the OME model of Ref. for the weak nonmesonic decay of hypernu- 
clei to the light of the new Nijmegen baryon-baryon potentials |^|. These strong interaction 



models influence the weak decay mechanism, not only through the coupling constants and 
form factors at the strong vertex involved in the AA^ NN reaction, but also via the PC 
piece of the weak vertex, obtained from a pole model, as well as from the corresponding 
correlated wave functions for the initial AA^ and flnal A^A^ states. 

We have paid an especial attention to disentangle the influence of the various model 
ingredients involved in the decay mechanism. In particular, we have found that the uncer- 
tainties tied to the strong coupling constants are a source of uncertainty on the nonmesonic 
weak decay observables. Employing different strong models of the Nijmegen group, and 
working consistently within each one, we estimate the decay rate of ^He to be in the range 
0.32 — 0.43 Fa, the neutron-to-proton ratio in the range 0.34 — 0.46 and the hypernuclear 
asymmetry parameter Ap around —0.68. In the case of ]^C, the decay rate turns out to 
be in the range 0.55 — 0.73 Fa, the neutron-to-proton ratio in the range 0.29 — 0.34 and the 
asymmetry Ap around 0.36. The new results for the nonmesonic rates compare favourably 
with the present experimental data. The ratio F„/Fp has increased with respect to our previ- 
ous works and it now lies practically within the lower side of the error band. The asymmetry 



for is also compatible with experiment ||3^ but that for ^He disagrees strongly from the 



recent experimental observation The latter work flnds a small and positive value for 
the elementary asymmetry parameter a a in AHe, while that for j^C is large and negative. 
Our meson-exchange model does not explain the present experimental differences and un- 
derstanding this issue is one of the current challenges, both experimental and theoretical, in 
the study of the weak decay of hypernuclei. 

We have found a tremendous influence on the weak decay observables from the way FSI 
are considered, especially in the case of total and partial decay rates. A phenomenological 
implementation of FSI effects, or not including them at all, gives rise to decay rates that 
differ by more than a factor of two, and to a neutron-to-proton ratio about 20% larger from 
what is obtained with the more realistic calculation that uses the proper A^A^ scattering 
wave function. The differences observed in the decay rates and the F„/Fp ratio are much 
larger than the uncertainties tied to the different strong interaction models commented 
above. Therefore, accurate calculations of the nonmesonic weak decay of hypernuclei demand 
a proper treatment of FSI effects through the solution of a T-matrix using realistic NN 
interactions. 
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APPENDIX 

The wave function describing the relative motion of two particles moving under the 
influence of a two-body potential V is obtained from the Lippmann-Schwinger equation: 

where | $) =| A;) is the free relative motion state, solution of the free Schrodinger equation 
with energy E , Hq \ ^) = E \ ^) , with Hq being the relative kinetic energy operator. The 
positive (negative) sign in front of the infinitesimal quantity irj arises from the requirement 
that the relative motion is free in the remote past (future). Alternatively, one can write 
the Lippmann-Schwinger equation in terms of the transition matrix T, defined such that 
T I $) = 1/ I ^W) ( (^(-) I -1/ = ($ I T), with the result 

I +(* I T^-^. (A.4) 

where the T operator obeys: 

Let's take for instance (representing the two nucleons in the final state) and project 
Eq. (|A.4| ) into coordinate space. Inserting a complete set of states in the r.h.s. of Eq. ( [A.4| ) 



and considering also the spin quantum numbers via two-particle coupled spin states, we find: 
{¥^\SMs I f) = {k,SMs I r) 

s^isJ E{k)-E{k>) + irj 

^ g-ifcr 

where (k, S Ms I r ) stands for the adjoint of the free plane wave, - — -irrr { SMs |. We per- 

(27r)'^/^ 

form a partial wave decomposition of the wave functions (\E'^ SMs \ r) and (k, S Ms \ r), 
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working in the coupled {LS)J representation, and, similarly to Eq. (18.119b) of Ref. 
we obtain: 



{¥f\SMs\f) 



E ( 

LL'S'JM 



^L'S',LS 



{k,r) 



X 



Ml 



where <Jijs'j{f') is the adjoint of the generalized spherical harmonic given by 



J2 {JM I L'M'LS'M's)YL>M',if) I S'M's) . 



(A.7) 



(A.8) 



The partial wave decomposition for the free wave function is simply obtained by replacing 
in Eq. (^) 



'^L's*,Lsik,r) ^ jL{kr)SLuSss' , 

where jL{kr) is the spherical Bessel function. 

A similar decomposition may be written down for the T-matrix elements, 

{k,SMs I T I k', S'M's) = E E E YLMAk)Yl,M',{k') 

JM LMl L'M'^ 

X {LMlSMs I JM){L'M'^S'M's \ JM) 
X {k{LS)JM I T I k'{L'S')JM) . 



(A.9) 



(A.IO) 



Inserting Eqs. (|A.7| ), and ( |A.10| ) into Eq. (|A.6| ), using Eqs. ( |A.8| ) and ( |A.9| ), and carrying out 
the angular integration, one obtains the following equation for the partial wave components 
of the correlated wave function 



'^L's*,Lsik,r) = jL{kr)6LL'Sss' 



+ / k'^dk 



,^,,j{k{LS)JM I T I k'{L'S')JM)jL'ik'r) 



E{k) - E{k') + ir] 
where the partial wave T-matrix elements fulfil the integral equation: 

{k{LS)JM I T I k'{L'S')JM) = {k{LS)JM \ V \ k'{L'S')JM) 
+ Yl I k"''dk 

S"L" 



(A.ll) 



„^^^^„{k{LS)JM\V\k"{L"S")JM) 



E{k) - E{k") + ir] 
X {k"{L"S")JM I T I k'{L'S')JM) . 



(A.12) 



Both equations ([A.ll|) and ( [A.12| ) are solved in momentum space following a numerical 
matrix inversion method described in Ref. [0. Isospin is easy to incorporate at the final 
step of our calculation by multiplying the resulting correlated wave function by the isospinor 
Xt"^ , allowing only for those cases that fulfil the antisymmetry requirement L + S + T = odd. 

An alternative but approximate way of implementing the A^A^ final state interactions is 
through the i^-matrix operator, which obeys: 
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K = V + VV 



E — Hn 



K 



(A.13) 



where V indicates the Cauchy principal value. Standing waves can then be obtained by 
following the same procedure of Eqs. ([A.6|) through (|A.12| ) but replacing the T-matrix 
elements by those of the real i^'-matrix. 

Finally, a simple way to include the effects of the strong interaction between the particles 
is to use a phenomenological correlation function, /phenl'"), such that 



fphen{r)jL{kr)6LL'SsS' 



(A.14) 



This has been the approach followed by several authors, and functions such as a Gaussian 
or of the type /phen('") = 1 — joiQc^) with = 3.93fm"^ have been used . 
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TABLES 



TT 


V 


K 


P 


UJ 


K* 


1750 


1750 


1789 


1232 


1310 


1649 



TABLE I. Cut-off values in MeV used in the present calculation for a FF of the type 
A^/(A^ + g^), which matches the exponential-type FF used in Ref. |28|. 



iHe 


r 

^ nm 


Lfi/Fp 




Ap 


TT 


0.438 


0.104 


0.397 


-0.282 


TT + K 


0.321 


0.286 


0.249 


-0.484 


all 


0.496 


0.226 


0.405 


-0.447 


A 


r 

-L nm 


Fn/Fp 




Ap 


TT 


0.771 


0.093 


0.705 


0.205 


TT + K 


0.558 


0.210 


0.461 


0.305 


all 


0.834 


0.181 


0.706 


0.275 



TABLE II. Weak decay observables for ^He and j^C including vr, vr + AT and all meson 
{7T + r] + K + p + uj + K*) contributions. The total and partial non-mesonic decay rates are 



in units of Fa = 3.8 x 10^ s~^. The strong NSC93 coupling constants |27| and Jiilich B cut-offs 
p6| in the (monopole) FF have been used. For the final NN wave function we used the solution of 
a T-matrix equation -Eq. (|A.6| )- with the NSC93 potential model. 
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iHe 


r 

a f 


a f 


a f 


a f 


vr 

■K + K 


0.424 0.425 
0.272 0.235 


0.086 0.086 
0.288 0.498 


0.390 0.391 
0.211 0.157 


-0.252 -0.252 
-0.572 -0.606 


all 


0.425 0.317 


0.343 0.457 


0.317 0.218 


-0.675 -0.682 


EXP: 


0.41 ±0.14 


181 


0.93 ± 0.55 [|1|] 


0.21 ±0.07 |18| 


0.24 ±0.22 |3€l 


12p 


r 

a f 


a f 


a f 


Ap 

a f 


vr 

■K + K 


0.762 0.751 
0.485 0.413 


0.078 0.079 
0.205 0.343 


0.707 0.696 
0.403 0.308 


0.169 0.171 
0.313 0.320 


all 


0.726 0.554 


0.288 0.341 


0.564 0.413 


0.358 0.367 


EXP: 


1.14 ±0.08 
0.89 ± 0.15 ±0. 
1.14 ±0.2 1 


37| 

03 [19| 
18| 


l-33i^i? |181 
1.87±0.59t;:32 [||] 

0.70 ±0.3 [20| 
0.52 ±0.16 


o.3ii[!i? m 


0.05 ± 0.530 



TABLE III. Weak decay observables for ^He and including vr, vr ± and all meson 
{'K + rj + K + p + ijO + K*) contributions. The total and partial non-mesonic decay rates are 
in units of T\ = 3.8 x 10^s~^. The strong NSC97a (left column) and NSC97f (right column) 
coupling constants and cut-offs in the FF (see Table |) have been used j^^. For the final NN wave 



function we used the solution of a T-matrix equation -Eq. ( A. 6 )- with the corresponding NSC97a 
and NSC97f potential models. 



^This number has been obtained dividing the experimental asymmetry, A = —0.01 ±0.10 |38], by 
a polarization of Py = —0.19 [39|. 
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iHe 


r 

nm 






Ap 


NSC97a 


0.320 


0.459 


0.219 


-0.680 


NSC97f 


0.317 


0.457 


0.218 


-0.682 


NSC93 


0.405 


0.483 


0.273 


-0.568 


Bonn B 


0.398 


0.484 


0.268 


-0.602 



TABLE IV. Weak decay observables for iHe including all meson {7r + r] + K + p + uj + K*) 
contributions. The total and partial non-mesonic decay rates are in units of Fa = 3.8 x 10^s~^. 
The strong NSC97f coupling constants and cut-offs in the FF (see Table 



have been used |28 |. 

For the final NN wave function we solve a T-matrix equation using different interaction models, as 
quoted in the table. 



iHe 


r 

nm 






Ap 


T 


0.317 


0.457 


0.218 


-0.682 


Re(T) 


0.490 


0.512 


0.324 


-0.655 


K 


0.475 


0.471 


0.323 


-0.650 


fphenir) 


0.766 


0.619 


0.473 


-0.671 


no FSI 


0.721 


0.614 


0.447 


-0.654 



TABLE V. Different approaches to the A^A^ final state wave function in the weak decay of iHe. 
The total and partial non-mesonic decay rates are in units of Fa = 3.8 x 10^s~^. The strong BB 
NSC97f model has been used. 
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FIGURES 




-0.2 



-0.4 



2 3 1 2 

r [fm] r [fm] 

FIG. 1. NN wave function in various channels for a relative momentum = 386 MeV/c, 
obtained from different models of the NN interaction: NSC97f (Hi (sohd line), NSC93 [|o| (dotted 
line) and Bonn B (dashed line). To avoid line overcrowding, only the real part of the wave 
function is shown in the plots. 
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-0.2 



-0.4 



2 3 1 2 

r [fm] r [fm] 

FIG. 2. NN wave function in various channels for a relative momentum pj. = 386 MeV/c, 
obtained from different treatments in the solution of the scattering equation: T-matrix (solid lines) , 
Real part of T-matrix (dotted line), K-matrix (dashed line) and phenomenological (dot-dashed 
line) . 
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